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We study non-Gaussian signatures on the cosmic microwave background (CMB) radiation pre- 
dicted within inflationary models with non-vacuum initial states for cosmological perturbations. 
The model incorporates a privileged scale, which implies the existence of a feature in the primordial 
power spectrum. This broken-scale-invariant model predicts a vanishing three-point correlation func- 
tion for the CMB temperature anisotropics (or any other odd-numbered-point correlation function) 
whilst an intrinsic non-Gaussian signature arises for any even-numbered-point correlation function. 
We thus focus on the first non- vanishing moment, the CMB four-point function at zero lag, namely 
the kurtosis, and compute its expected value for different locations of the primordial feature in the 
spectrum, as suggested in the literature to conform to observations of large scale structure. The 
excess kurtosis is found to be negative and the signal to noise ratio for the dimensionless excess 
kurtosis parameter is equal to \S/N\ ~ 4 x 10~ 4 , almost independently of the free parameters of 
the model. This signature turns out to be undetectable. We conclude that, subject to current tests, 
Gaussianity is a generic property of single field inflationary models. The only uncertainty concern- 
ing this prediction is that the effect of back-reaction has not yet been properly incorporated. The 
implications for the trans-Planckian problem of inflation are also briefly discussed. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

The theory of inflation is presently the most appealing candidate for describing the early universe. Inflation 
essentially consists of a phase of accelerated expansion which took place at a very high energy scale. One of the main 
reasons for such an appeal is the fact that inflation is deeply rooted in the basic principles of general relativity and 
field theory, which are well-tested theories. It is because all the forms of energy gravitate in general relativity that 
one of them, the pressure, which can be negative in field theory, is able to cause the acceleration in the expansion 
of the universe. In addition, when the principles of quantum mechanics are taken into account, inflation provides 
a natural explanation for the origin of the large scale structures and the associated temperature anisotropies in the 
Cosmic Microwave Background (CMB) radiation [1]. 

Inflation makes four key predictions: (i) the curvature of the space-like sections vanishes, i.e. the total energy density, 
relative to the critical density, is f^o = 1, (ii) the power spectrum of density fluctuations is almost scale invariant, 
i.e. its spectral index is n s ~ 1, (iii) there is a background of primordial gravitational waves (which is also scale 
invariant), and (iv) the statistical properties of the CMB are Gaussian. In this article we focus on the last prediction 
and investigate whether it is a robust and generic property of inflationary models. The statistical properties of the 
CMB will be measured with high accuracy by the MAP and Planck satellites [2] . So far the preliminary measurements 
of the three- and four-point correlation functions [3, 4, 5] seem to be consistent with Gaussianity. 

The fact that the statistical properties of the CMB are Gaussian can be directly traced back to the common 
assumption that the quantum fluctuations of the inflaton field are placed in the vacuum state [6]. Therefore, in 
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order to answer the above question, one has to investigate which kind of non-Gaussianity shows up if the vacuum 
state assumption is relaxed [7, 8, 9, 10]. In particular, one crucial point is to study whether this modification yields 
a detectable signal for future CMB or large-scale structure observations. Let us also notice that there exist other 
mechanisms to produce non-Gaussianity within the framework of inflation. Some of them have been studied in 
Rcfs. [11, 12]. 

Assuming that the quantum state of the perturbations is a non-vacuum state immediately leads to the following 
difficulty: non-vacuum initial states imply, in general, a large energy density of inflaton field quanta, not of a cosmo- 
logical term type [13]. In other words, generically, if the initial state is not the vacuum then there is a back-reaction 
problem that could upset the inflationary phase. However, as we will argue below, one cannot directly conclude that 
this would prevent inflation from occurring altogether because, without a detailed calculation, it is difficult to guess 
what the back-reaction effect on the background would be. Such a detailed calculation is in principle possible by 
means of the formalism developed in Ref. [14]. To our knowledge, such a computation has never been performed. The 
calculation of second order effects is clearly a complicated issue and is still the subject of discussions in the literature, 
see [15] for example. Moreover there exist situations where it can be avoided, and this is in fact the case if the number 
of e-folds is not too large. In this article, we will not address the general question mentioned above but will rather 
concentrate on the more modest aim of calculating the non-Gaussianity in a case where the back-reaction problem is 
not too severe, hoping in this way to capture some features of the real situation. 

There exist other arguments to study the non-Gaussianity that arises from a non-vacuum state. One of these is 
the so-called trans-Planckian problem of inflation [16]: the quantum fluctuations are typically generated from sub- 
Planckian scales and therefore the predictions of inflation depend in fact on hidden assumptions about the physics 
on length scales smaller than the Planck scale. However, it has recently been shown that inflation is robust to some 
changes of the standard laws of physics beyond the Planck scale. More precisely, inflation is robust to a modification of 
the dispersion relation, at least if those changes are not too drastic, in practice if the Wigner-Kramer-Brillouin (WKB) 
evolution of the cosmological perturbations is preserved. However, modeling trans-Planckian physics by a change in 
the dispersion relation is clearly ad-hoc. Therefore, it is interesting to consider other possibilities; for example, one 
could imagine that the inflaton field emerges from the trans-Planckian regime in a non- vacuum state. Non-Gaussianity 
would then be, in this case, a signature of non standard physics and it seems to us interesting to quantify this effect. 
Let us note that similar ideas have been suggested in Ref. [17] in a slightly different context. Let us also remark that 
it has been shown recently in Ref. [18] that placing the cosmological perturbations in a non-vacuum state would lead 
to possible observable effects, for instance a modification of the consistency check of inflation. 

Another motivation for calculating non-Gaussianity when the initial state is not the vacuum is that this model 
could be used to test the methods that are being developed to detect non-Gaussianity in the future CMB maps. 
There have been approaches based on the nth order moments or the cumulants of the temperature distribution [19], 
the n-point correlation functions or their spherical harmonic transforms [20] , and also works based on the detection of 
gradients in the wavelet space [21], to mention just a few methods. In the last approach, namely the wavelet analysis 
of a signal [21], the test maps which were employed were characterized by a non-skewed non-Gaussian distribution. 
Therefore, any non-Gaussianity was indicated by a non-zero excess kurtosis of the coefficients associated with the 
gradients of the signal. As we will show, this is exactly our case. This wavelet analysis of a signal was then applied [22] 
to search for the CMB non-Gaussian signatures. More precisely, these authors investigated the detectability of a non- 
Gaussian signal induced by secondary anisotropies, while assuming Gaussian-distributed primary anisotropies. Their 
method [22] is unable to detect such non-Gaussianity for the MAP-like instrumental configuration while it can do it 
for Planck-like capabilities. 

From the theoretical point of view the simplest way to generalize the vacuum initial state, which contains no 
privileged scale, is to consider an initial state with a built-in characteristic scale, kb [8]. Here, we will consider a non- 
vacuum state which is simpler and more generic than the one considered in this previous work. Several observables 
can be used to constrain the parameter space. A first possibility is to use the CMB anisotropy multipolc moments 
to constrain the number of quanta n around the privileged scale. It has been shown in Ref. [8] that, typically, this 
number cannot be large and in the present article we will always consider that n is a few. With the recent release of 
the BOOMERanG [23], MAXIMA [24] and DASI [25] data, which revealed the existence of a first acoustic peak in 
the angular power spectrum at I <~ 200, followed by a second acoustic peak located at I ~ 500 and an evidence for a 
third peak, one can hope to obtain stronger constraints on kb and n very soon. Of course, another observable which 
can also be used is the matter density power spectrum. We will compare the predictions of our model for different 
cosmologies with the result of recent observations below. 

The model we are studying here belongs to a class with a Broken Scale Invariant (BSI) power spectrum for the 
matter density. Such a primordial spectrum could also be generated [26] during an inflationary era where the inflaton 
potential is endowed with steps, e.g., induced by a spontaneous symmetry breaking phase transition. The main 
motivation behind this class of models comes from Abell/ACO galaxy cluster redshift surveys which indicate [27] that 
the matter power spectrum seems to contain large amplitude features close to a scale of 100/i _1 Mpc (see however 
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[28]). In support of this finding arc the preliminary results of the recently released [29, 30] power spectrum analysis of 
the rcdshift surveys of Quasi-Stellar Objects (QSOs). Using the 10k catalogue from the 2dF QSO Rcdshift Survey [31], 
it has been tentatively identified [29] a "spike" feature at a scale w 90/i _1 Mpc (rs 65/i _1 Mpc) assuming a Lambda 
contribution S}\ = 0.7 and an ordinary matter contribution £l m = 0.3 (respectively, = and f2 m = 1.0). Provided 
this feature is confirmed, it might also have originated from acoustic oscillations in the tightly-coupled baryon-radiation 
fluid prior to decoupling. Using the CMBFAST code it was found [29] that this spike in the spectrum is seen at a 
> 25 per cent smaller wavenumber than the second acoustic peak, while higher values of the baryon contribution fib 
may be needed to fit the amplitude of this feature. If we interpret this feature as originating from the primordial 
spectrum then, in order to be consistent with observations, the preferred scale fcb must lie way below the horizon 
today, possibly at a scale corresponding to the turn of the power spectrum [27] or at the scale matching the first 
acoustic peak of the CMB temperature anisotropics [32]. One then sees that the presently available data already 
restricts the parameter space for the quantities fcb and n [33]. 

For the class of models which contain a preferred scale, a generic prediction is that the three-point correlation 
function vanishes (as well as any higher-order odd-point function), whereas the four-point (as well as any higher-order 
even-point) correlation function no longer satisfies the relation 



which is typical of Gaussian statistics. Since the third-order moment (the skewness) vanishes, a first step is to 
calculate the fourth-order statistics (the kurtosis). It is interesting to perform this calculation for very large (COBE- 
size) angular scales, for which one can be confident that the source of non-Gaussianity is primordial. On the other 
hand, if the non-Gaussian signature was calculated on intermediate scales, a stronger signal would be obtained; 
however, in that case the secondary sources would be more difficult to subtract and thus the transparency of the effect 
would be compromised. To quantify the relevant amplitude of the signal, the excess kurtosis should be compared 
with its cosmic variance. This was computed, e.g., in Ref. [34], for a Gaussian field. Although, strictly speaking, one 
should compute the cosmic variance for the actual case and not rely on a mildly non-Gaussian analysis, the actual 
smallness of the obtained signal largely justifies our approach. 

We organize the rest of the paper as follows: in Section II, we discuss in detail the argument developed in Ref. [13] 
regarding the back-reaction problem. We show that any theory with a non-vacuum initial state has to face this issue. 
However, we also argue that it is not clear at all whether inflation will be prevented in this context. In Section III, 
we discuss our choice of a non-vacuum initial state for cosmological perturbations of quantum-mechanical origin and 
we give some basic formulas for the two-point function. We calculate the CMB angular correlation function and the 
associated matter power-spectrum for our choice of non-vacuum initial states, comparing the latter against current 
observations. In Section IV, we calculate the angular four-point correlation function and the related CMB excess 
kurtosis, while in Section V we discuss our results explicitly and present both analytical and full numerical estimates 
of the normalized excess kurtosis for a typical case. In this Section we also compare this non-Gaussian signal with 
its corresponding cosmic variance. We round up with our conclusions in Section VII. In this paper, we use units such 
that c = 1. 



II. INITIAL STATE FOR THE COSMOLOGICAL PERTURBATIONS AND THE BACK-REACTION 

PROBLEM 



In this Section we discuss the relevance of non-vacuum initial states for cosmological quantum perturbations. The 
argument of Ref. [13] is based on the calculation of the energy density of the perturbed inflaton scalar field in a given 
non-vacuum initial state. Since the perturbed inflaton and the Bardcen potential are linked through the Einstein 
equations, it is clear that they should be placed in the same quantum state. Let us consider a quantum scalar field 
living in a (spatially flat) Friedmann-Lemaitre-Roberston- Walker background. The expression of the corresponding 
operator reads 



^(77, x) = 
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where Ck(?7i) and c k (r]i) are the annihilation and creation operators (respectively) satisfying the commutation relation 
[ck, 4>] = 5(k— p), and where a(r]) is the scale factor depending on conformal time 77. The equation of motion for the 
mode function Hk(v) can be written as [35, 36, 37] 

Mfc+(fc 2 -^)^ = 0, (3) 
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where "primes" stand for derivatives with respect to conformal time. The above is the characteristic equation of 
a parametric oscillator whose time-dependent frequency depends on the scale factor and its derivative. The energy 
density and pressure for a scalar field are given by the following expressions 
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Let us now calculate the energy and pressure in a state characterized by a distribution n(k) (giving the number n of 
quanta with comoving wave-number k) for a free (i.e. V = 0) field. Let us denote such a state by \n(k)). Using some 
simple algebra it is easy to find 
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We can evaluate these quantities in the high-frequency regime and take p,k ~ exp[— ik(r\ — rfi)], where r)i is some given 
initial conformal time. We get 
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Several comments are in order at this point. Firstly, the lower limit of the integral is certainly not zero because at 
some fixed time, k — > corresponds to modes outside the horizon. So if we evaluate the previous integral at time 
r] then we should only integrate over those modes whose wavelength is smaller than the Hubble radius. But in the 
infra-red sector, the integral is finite and so the contributions of those modes will be small. Therefore, in practice 
we can keep a vanishing lower bound. Secondly, the first term of each expression is the contribution of the vacuum, 
i.e., is present even if n(k) = 0. This is clearly divergent in the ultra-violet regime. At this point, one should adopt 
a regularization procedure (in curved space-time). Once this infinite vacuum contribution is subtracted out, our 
renormalized expressions for the density and pressure in the \n(k)) state read 
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For a well-behaved distribution function n(k) this result is finite. Thirdly, the perturbed inflaton (scalar) particles 
behave as radiation, as clearly indicated by the equation of state p = (l/3)p and as could have been guessed from the 
beginning since the scalar field studied is free. To go further, we need to specify the function n(k). If we assume that 
the distribution n(k) is peaked around a value fcb, it can be approximated by a constant distribution of n quanta, 
with n(fcb) — n, in the interval [fcb — Afc, fcb + Afc] centered around fcb. If the interval is not too large, i.e. Afc -c fcb 
then, at first order in Afc/fcb, we get 
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where iV e is the number of e-folds counted back from the time of exit, see Fig 1. The time of exit is determined by 
the condition fc p h ys = k/a ~ i?; n f, where H ln f is the Hubble parameter during inflation. It is simply related to the 
scale of inflation, M; n f, by the relation £f; n f ~ M 2 n{ /m pl . We have also assumed that, during inflation, the scale factor 
behaves as a(t) oc exp(iJ; n ft). From Eq. (12), we see that the back-reaction problem occurs when one goes back in 
time since the energy density of the quanta scales as ~ 1/a 4 . In this case, the number of e-folds N e increases and the 
quantity (n(k)\p\n(k)} raises. This calculation is valid as long as (n(k)\p\n(k)) < p- m i = m^H 2 ^. When these two 
quantities are equal, the energy density of the fluctuations is equal to the energy density of the background and the 
linear theory breaks down. This happens for 7V C = N hr such that 
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where we have assumed nA£;/(7r 2 fcb) — 0(1). Interestingly enough, this number does not depend on the scale k but 
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FIG. 1: Sketch of the evolution of the physical size in an inflationary universe where the cosmological perturbations are placed 
in a non- vacuum state characterized by the distribution n(k). 



only on the Hubble radius during inflation H m f. This means that for each scale considered separately, the back- 
reaction problem starts to be important after the same number of e-folds A c = A br counted back from horizon exit 
(this is why in Fig. 1, one has A br = A br for the two different scales Ai and A3). If, for instance, we consider the 
case where inflation takes place at GUT scales, M- m f ~ 10 16 GeV, then iJi n f — 10 13 GeV and one obtains A br ~ 7 in 
agreement with the estimates of Ref. [13]. If the distribution n(k) is not strongly peaked around a particular scale but 
is rather spread over a large interval, it is clear that the important mode of the problem is, very roughly speaking, the 
populated smallest scale (i.e. ~ Ai in Fig. 1). In the following, this scale is denoted by A pop . The value of this scale 
clearly depends upon the form of the distribution n(k). As it can be seen in Fig. 1, A pop is the scale for which the 
back-reaction problem shows up first, as we go backward in time, since the other modes with larger wavelengths have 
not yet penetrated deeply into the horizon and therefore do not yet face a back-reaction problem. As a consequence, 
this scale determines the total number of e-folds during inflation without a back-reaction problem. 

Once the number A br has been calculated, the total number of e-folds during inflation without a back-reaction 
problem is a priori fixed. It remains to be checked whether this number is still sufficient to solve the usual problems of 
the hot big bang model. We now turn to this question. Let N» (A) be the number of e-folds, for a given scale A, between 
horizon exit during inflation and the beginning of the radiation era, see Fig. 1. The total number of e-folds of inflation 
without a back-reaction problem is then A^ n f = A br + A*. The number A* is given by A* (A) = m(ao/a*) — A r — A m , 
where ao and a* (A) are the scale factor at present time and at first horizon crossing, respectively. The quantities A r 
and A m are the number of e-folds during the radiation and matter dominated epochs. The ratio ao/a* is given by 
{^/^n)Hi n {/H , where l H is the present day Hubble radius and H is the present value of the Hubble parameter given 
by H /m pl ~ 10~ 61 . The quantities A r and A m are given by A r = ln(T r h/T e q) and A m ~ ln(z cq ) ~ 9. T r h is the 
reheating temperature which can be expressed as T r h ~ (Tm^) 1 ^ 2 where T is the decay width of the inflaton. For 
consistency, one must have M; n f > T r h- T eq is the temperature at equivalence between radiation and matter and its 
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The quantity A* (A) can be expressed as 
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From now on, in order to simplify the discussion, we assume that the decay width of the inflaton field is such that 
T r h — Mi„f. Under these conditions, the usual problems are solved if the number Aj n f is such that Ai n f(A pop ) ~ 
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ln(A pop /^ H ) + 29 x In 10 > —4 + lnz cn( j, where the quantity z cn d is the redshift at which the standard evolution (hot 
big bang model) starts. It is linked to the reheating temperature by the relation log 10 (z cn( j) — 32 + log 10 (T r i 1 /m pl ). 
This gives a constraint on the scale of inflation, namely 

logl0 (^i) <21ogl0 (W)_, 5 . (15) 

It is known that inflation can take place between the TeV scale and the Planck scale which amounts to —32 < 
log 10 (.ffinf //rip, ) < 0. We see that the constraint given by Eq. (15) is not too restrictive. In particular, if we take 
Apop = 0.1£ H and Hi n f = 10 13 GeV, it is satisfied. However, if we decrease the scale A pop , the constraint becomes more 
restrictive. The constraint derived in the present article appears to be less restrictive than in Ref. [13] because we do 
not assume that all scales are populated. 

Another condition must be taken into account. We have seen that the duration of inflation without a back- reaction 
problem is determined by the evolution of A pop . However, at the time at which the back-reaction problem shows 
up, one must also check that all the scales of astrophysical interest today were inside the horizon so that physically 
meaningful initial conditions can be chosen. This property is one of the most important advantages of the inflationary 
scenario. If we say that the largest scale of interest today is the horizon, this condition is equivalent to 

JV.(4) < JV„(A pop ) + N hl N h * > ]n(^-) . (16) 

V A pop / 

This condition is also not very restrictive, especially for large scales. As previously, the condition can be more 
restrictive of one wants to populate smaller scales. 

Before concluding this section, a last comment is in order. What is actually shown above is that, roughly N hr 
c-folds before the relevant mode left the horizon, we face a back-reaction problem, as the energy density of the 
perturbation (n(k)\p\n(k)} becomes of the same order of magnitude as the background p. So, before concluding that 
non- vacuum initial states may or may not turn off the inflationary phase, one should calculate the back-reaction effect, 
i.e., extend the present framework to second order as it was done in Ref. [14]. To our knowledge, this analysis is still 
to be performed. Moreover, even if we take the most pessimistic position, that is, one in which we assume that the 
back-reaction of the perturbations on the background energy density prevents the inflationary phase, there still exist 
models of inflation where the previous difficulties do not show up. Therefore, in the most pessimistic situation, there 
is still a hope to reconcile non-vacuum initial states with inflation. Admittedly, the price to pay is a fine-tuning of 
the free parameters describing inflation and/or the non- vacuum state. 

III. TWO-POINT CORRELATION FUNCTION FOR NON- VACUUM INITIAL STATES 

A. General expressions 

We now turn to consider the non-vacuum states for the cosmological perturbations of quantum mechanical origin. 
Let V{a) be a domain in momentum space, such that if k is between and a, the domain T>{a) is filled by n quanta, 
while otherwise D contains nothing. Let us note that this domain is slightly different from the one considered in 
Ref. [8]. The state |*i(c7,n)> is defined by 

(J )n 

|*i(a,n))= I] 4=-|°0 ® |0p>= Ink) (g) |0 P >. (17) 

keX>(<x) p£X>(cr) k£X>(cr) p£X>(<r) 

The state |nk) is an n-particle state satisfying, at conformal time r] = 771 : Ck|nk) = y/n\(n — l)k) and cjjnk) = 
y/n + l\(n + l)k)- We have the following property 1 

{9 1 {*,n)\* 1 (o',ri))=5((T-o')8 nn ,. (18) 

It is clear from the definition of the state that the transition between the empty and the filled modes is sharp. 
In order to "smooth out" the state we consider a state ^2} as a quantum superposition of In doing so, 



This normalization is in agreement with Eq. (2.25) of Ref. [38]. 
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we introduce an, a priori, arbitrary function g(a; fcb) of a. The definition of the state \^2{n, fcb)) is 

r+oo 

|* 2 (n,fcb))=/ do-</(tr;fcb)|*i(tr,n)), (19) 
Jo 

where g(a; fcb) is a given function which defines the privileged scale fcb- We assume that the state is normalized and 
therefore Jq + °° g 2 (o~', fcb)der = 1. In the state |5'i(cr, n)}, for any domain T> one has [8]: 

(tfiforOlcpCql^ifon)) = (*i(tr,n)|44|*i(o-,n))=0, (20) 
(*i(tr,n)|cp4|*i(ff,n)) = n<5(q e ©)5(p - q) + 5(p - q), (21) 
(* 1 (<7,n)| C t,c q |'I' 1 (<7,n)) = nJ(q € ©)<y(p - q). (22) 

In these formulas, <5(q € 2?) is a function that is equal to 1 if q £ V and otherwise. These relations will be employed 
in the sequel for the computation of the CMB temperature anisotropics for the different non-vacuum initial states. 



B. Two-point correlation function of the CMB temperature anisotropy 



The spherical harmonic expansion of the cosmic microwave background temperature anisotropy, as a function of 
angular position, is given by 
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The Wi stands for the ^-dependent window function of the particular experiment. In the work presented here, we are 
interested in a non-Gaussian signature of primordial origin. We are thus focusing on large angular scales, for which 
the main contribution to the temperature anisotropy is given by the Sachs- Wolfe effect, implying 



^r(e) - ^$[?7is S ,e(77o - r?i ss )], 



(24) 



where $(77, x) is the Bardeen potential, while 770 and 7?i ss denote respectively the conformal times now and at the last 
scattering surface. Note that the previous expression is only valid for the standard Cold Dark Matter model (sCDM). 
In the following, we will also be interested in the case where a cosmological constant is present (ACDM model) since 
this seems to be favored by recent observations. Then, the integrated Sachs- Wolfe effect plays a non- negligible role 
on large scales and the expression giving the temperature fluctuations is not as simple as the previous one. 

In the theory of cosmological perturbations of quantum mechanical origin, the Bardeen variable becomes an oper- 
ator, and its expression can be written as [8] 



$(r/,x) = 
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where £pi = (Gft) 1 / 2 is the Planck length. In the following, wc will consider the class of models of power- law 
inflation since the power spectrum of the fluctuations is then explicitly known. In this case, the scale factor reads 
a(v) = 4M 1+/3 , where (3 < -2 is a priori a free parameter. However, in order to obtain an almost scale-invariant 
spectrum, (3 should be close to —2. In the previous expression of the scale factor, the quantity £q has the dimension 
of a length and is equal to the Hubble radius during inflation if (3 = —2. The parameter £ also appears in Eq. (25). 
The factor 3/(4tt) in that equation has been introduced for future convenience: the factor 3 will cancel the 1/3 in 
the Sachs- Wolfe formula and the factor l/(4ir) will cancel the factor 4tt appearing when the complex exponentials are 
expressed in terms of Bessel functions and spherical harmonics. The mode function fk(f]) °f the Bardeen operator is 
related to the mode function /J-k(v) °f t nc perturbed inflaton through the perturbed Einstein equations. In the case 
of power-law inflation and in the long wavelength limit, the function fk(rj) is given in terms of the amplitude A s and 
the spectral index n s of the induced density perturbations by 



k 3 \fkf 



k 



Using the Rayleigh equation and the completeness relation for the spherical harmonics we have 



exp 



i(Vo -f/i ss )k-e 



= 4tt JVj/[fc(»») - il\ss)}Y; m (k)Y lm (e) 



(26) 



(27) 
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where jt denotes the spherical Bessel function of order t. Equations (23), (24), (25) and (27) imply 



Ipi m/2 
to 



dk 



Ck(m)fk(v) + c -k( J 7i)/fc( ? 7) 



je[k(vo - r? lss )]r/ m (k) . 



(28) 



At this point we need to somehow restrict the shape of the domain T>. We assume that the domain only restricts the 
modulus of the vectors, while it does not act on their direction. Then, from Eq. (28), one deduces 



I 2 f +oa dk 

(*i(o-,n)|a« imi a; am2 |*i(<r,n)) = 5 tlt3 8 mima ^ y & [k( m - ^ SB )]fc 3 |/ fc | 2 [l + 2n6(k e V)\ 



(2 
f 2 



with 



D { p{a) = I fi[k( m - m ss )]k'Vk\ 2 ^ = ^A s I J 2 +l/2 (k)k n ^dk ee ■-A.Df'ia) , 



C ll +2nD ( e 1 i ) (<i) 



,dk 



til2 u mim2 j 



k 



'1+1/2 



i.ns — 3. 



(29) 



(30) 



where Ji(z) is an ordinary Bessel function of order £. In the last equality and in what follows we take 770 — rjiss = 1- 
The amplitude A s and the spectral index n s are defined by Eq. (26). Thus, the multipole moments in the state 
l^i), are given by 



C { 1\(t) =C l + 2nDY\cT) 



(31) 



where Ct is the "standard" multipole, i.e., the multipole obtained in the case where the quantum state is the vacuum, 
i.e., n = 0. Let us calculate the same quantity in the state l^)- Performing a similar analysis as the above one, we 
find 



1*2(1, fcb)> 



£ 2 

l 

r>+oo 



r + °° c\k 
J jt[k{^-rh sa )]k s \h\ 2 T 

r-\-oo per 

+2n dag 2 (a;k h ) j\ [k{ m - r?i ss )]fc 3 |/ fc | 
Jo Jo 



.dk 

T 



(32) 



Defining g 2 (a; k\>) = dh/da [we will see below that this function h, actually /i(fcb), cannot be arbitrary] and integrating 
by parts leads to 



£2 r+00 f 

(^2(n,k h )\a eimi a* e2m2 \^ 2 (n, k h )} = 8 lll2 5 mim2 -^ J j\ [k(r) - ?7i ss )]fc 3 |/ fc | 2 |l + 2n/i(oo) 



hjk) 
h(oo) 



dk 
~k 



p2 
£P1 

l 



(2) 



C tl + 2nD\\ 



8e 1 e 2 8 mir 



(33) 



with 



D 



(2) 



h(oo) 



Ji[k(Vo - Vis 



h{k) 
h(oo) 



fc3|/fe|2 T = \ As J J?+i/2(k)h(k)k ns - 3 dk EE ^A S D {2) , (34) 



where h(k) = h(oo)[l — h(k)/h(oo)]. To perform this calculation, we have not assumed anything on h(oo) or h(0). We 
see that the relation g 2 (k) = dh/dk requires the function h(k) to be monotonically increasing with k. It is interesting 
that, already at this stage of the calculations, very stringent conditions are required on the function h(k) which is 
therefore not arbitrary. This implies that the function h(k) which appears in the correction to the multipole moments is 
always positive, vanishes at infinity and is monotonically decreasing with k. An explicit profile for h{k) is given below. 
In addition, the state \^2{n,kb)) must be normalized, which amounts to take, see Section III, g 2 (a;kb)da = 1. 
Using the definition of the function g 2 , we easily find 



ft(oo) - h(0) = 1 => h(0) = 1. 
The total power spectrum of the Bardeen potential can be written as 



/■■ :i |<I>/,| 2 x '{' ' 2>//)(-x) 



1 



hjk) 
h(oo) 



A s k n s- 1 [l + 2nh(k)}. 



(35) 



(36) 
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The exact proportionality coefficient is derived below. Observations indicate that n s ~ 1 and for simplicity we will 
take n s = 1. Then, if the function h(k) is chosen such that it contains a preferred scale, see the Introduction, and such 
that it is approximatively constant on both sides, the model becomes very similar to the one presented in Ref. [7] for, 
in the notation of that article, p > 1. In Ref. [7] the allowed range of parameters is 0.8 < p < 1.7 with an especially 
good agreement for the inverted step p < 1. In our case, another difference consists in the fact that the oscillations in 
the spectrum studied in Ref. [7] are not present due to the monotony condition on the function h(k). In the following, 
in order to perform concrete calculations, we will choose an analytical form for h(k) which mimics the behavior of the 
spectrum considered in Ref. [7] with p > 1. Interestingly enough, we will see that the final result does not strongly 
depend on the values of the free parameters that describe the function h(k). As we have seen previously, we can write 
the multipole moments in the state ^2) as 



Cf } =C t + 2nD<f } 



(37) 



Substituting the well-known expression for the Ci's and the definition of D^p given by Eq. (34), one finds that the 
coefficients C^p are given by 



1 r(3-n s )r[i + (n s -i)/2] 

2 1 2 3 -« S T 2 [(4 - n s )/2]T[£ - (n s - 5)/2] 



+ 2nD 



(2; 



(38) 



As a next step, one has to normalize the spectrum or, in other words, we need to determine the value of A$. We 

choose to use the value of <5 rms -F 
satellite. The quadrupole is then 



choose to use the value of Q rms _ PS = T [5Cf ) W 2 2 /(47r)] 1 / 2 (^ P1 /4) - 18fiK with T = 2.7K measured by the COBE 



a 



'-'> -a ~f 1 

3 2\2 3 -™s T 2 [(4-n s )/2]r[4- (n, 



1 "jr[2 + („, - i.)/2] +2n5 ( 2) | 



8 Qrms-PS ^0 



l)/2] 



5W 2 



T 2 



£ 2 



A- + 2nDf 

07T 



(39) 



where, in order to establish the last relation, and hereafter, we have assumed n s 
expressed in terms of band powers 5Tg defined by 



1. The measurements arc often 



8T e = 



'+1)4^(2) 



27T l\ ^ 



(40) 



For n s = 1, Eq. (39) can be simplified and the band powers defined by Eq. (40) lead to 



ST e 



Q 



rms-PS 



W 2 



1 + 2nir£(e + \)D 



(2) 



\ 1 + UnirD^ 



(2) 



(41) 



The n-dependence in the above expression is the correction due to the non- vacuum initial state. One can easily check 
that if n = the corresponding band powers are constant at large scales, a property which is well-known. 

Finally, we calculate the two-point correlation function at zero lag in the state 1^2) ■ Using Eqs. (23), (33), (37), 
the second moment, jU 2; 01 the distribution is given by 



M2 



jr(e) 



ll , ^2^+1 ^ (2) 



(I ^ 



47T 



(42) 



Once we have reached this point, an obvious first thing to do is to check that the two-point correlation function 
calculated above is consistent with present observations. 



C. Comparison with observations 

Among the available observations that one can use to check the predictions of theoretical models, two are key in 
cosmology: the CMB anisotropy and the matter-density power spectra. Since the initial spectrum is very similar to 
the one considered in Ref. [7], it is clear that the multipole moments and the matter power spectrum will also be 
similar to the ones obtained in that article. This already guarantees that there will be no clash with the observations. 
Therefore, we will not study in details all the predictions that can be done from the two-point correlation function 
since our main purpose in this article is to calculate the non-Gaussianity which is a clear specific signature of a non 



10 



vacuum state (in Ref. [7], the CMB statistics is Gaussian). Here, we just calculate the matter power spectrum to 
demonstrate that it fits reasonably well the available astrophysical observations for some values of the free parameters. 
In addition, this illustrates well the fact that, using the available observations, we can already put some constraints 
on the free parameters. In Ref. [8], although the model considered was slightly different, the multipole moments were 
computed and shown to be in agreement with the data if the number of quanta is a few. Therefore, having given all 
these reasons, it seems logical to concentrate in the present article on the matter power spectrum. 



1. Choice of the weight function 
We first need to choose the function h(k) such that it satisfies the conditions described above. A simple ansatz is 

h(k) = A + Btanhfaln-^- ] . (43) 



fc b 

In this equation, fcb represents the privileged (comoving) wave-number and a is a parameter which controls the 
sharpness of the function h(k). The argument of the hyperbolic tangent is expressed in terms of the logarithm of the 
wave- number in order to guarantee that fc € [0, +oo[, see, e.g., Ref. [32]. A and B are two coefficients that we are 
going to fix now. We have h(0) = A — B and h(oo) = A + B. Therefore, the requirement that the state be normalized 
translates into the condition B = 1/2. In fact, it is easy to see that h(k) does not depend on A. The function h(k) 
can be written as 



, ( , k 
1 — tanh am — 

V fcv 



(44) 



We can easily check that h(Q) = 1 and that h(oo) = 0. To make the connection with the literature, note that 
we are not taking ao = 1 today. Rather, ao — 2^h ~ 6000 h~ 1 Mpc which implies that a preferred scale located at 
0.004 h Mpc -1 [32] corresponds to a comoving fcb = 24 in our case, while a preferred scale located at 0.052 h Mpc 1 [27] 
yields fcb = 312 in our units. In order to illustrate the different forms of the initial spectrum, the function h(k) is 
represented in Fig. 2. 




FIG. 2: The function h(k) for a = 2 (continuous line), a = 5 (dotted line) and a = 10 (dashed line). 
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At this point, we need to re-investigate the back- reaction problem described before but now for the two non- vacuum 
states l^i) and | - In particular, we are going to see that we can put some constraints on the free parameters only 
from theoretical considerations. An analogous analysis to the one performed in Section II, now for the state 
leads to 



(*i|p|*i> 



1 



8vr 3 a 4 



S^a 4 I 



+°° dfc 
fc 

+°° dfc 
fc 



fc 2 [l + 2n<5(kel?)] 
fc 2 [l + 2n<5(kel?)] 



PfcPfc 



"(Pfc/4* + /4/4) + 



- 



PfcPfc 



(45) 
(46) 



Again, as in Section II, the term which is not proportional to n should be subtracted since it represents the vacuum 
contribution. Now, for the state \^2), one gets 



(*2|p|*2) 



00 /*+oo 



/' + OO 

L ' 



dada'g* (a) 9 (a')(*i (a, n)\p\^ l (a', n)) 



d<7| fl (a)| 2 <*i(<r, n)|/9|*i(^, ")}, 



where we used the fact that p does not act on a. In the high-frequency regime, one obtains 

n r°° 



(*2bl* 2 ) = (o|p|o) 



1 



3 167r 2 a 4 



^\g(a)\ 2 da. 



(47) 
(48) 

(49) 
(50) 



Clearly, the equation of state is still radiation. A comparison with Sec. II also shows that the distribution of quanta is 
now known explicitly and is given by n(fc) = n\g(k; fcb)| 2 fc/8. Using |g(c)| 2 = dh/da, the distribution function reads 



/i \ an 1 -2 
n(k) = — cosh 

v ; 16 



<\ In ( A 

fcb 



(51) 



This function is represented in Fig. 3 for different values of the parameter a. It is clearly peaked around the central 
value fc = fcb- Finally, we end up with the following energy contribution from our state (after having removed the 
vacuum contribution) 



(* 2 |p|*2> = 



167r 2 a 4 2 



1 — tanh a In 



fcb 



da 



fc 4 



16ir 2 a 4 2 



e 4v/a 

cosh 2 v 



dv, 



(52) 



where we have used the change of variable v = adn(er/fcb). Separating the integral into two pieces, it is easy to show 
that 



cosh(4w/o;) 
cosh 2 v 



dv = 



16^ 2 o 4 b 



k£B 1 + -, 1- 



a 



a 



167r 2 a 4 b 



ktr 1 + 



r 1- 



Q 



(53) 



where we have used Eqs. (3.512.1) and (8.384.1) of Ref. [39] and where B(x,y) = ^^{l - t) v ~ l dt is the Euler's 
integral of the first kind, see Eq. (8.380.1) of Ref. [39]. In the above equation T is the Euler's integral of the second 
kind. These expressions are well-defined if a > 2; it is interesting to see that we can obtain constraints on the free 
parameters of our model just from the requirement that the energy be finite. We see in Fig. 2 that it corresponds 
to a situation where the function is almost a step. Repeating the same reasoning as before, we obtain the following 
constraint: 



ffinf 



— In n — In 



r 1 + 



Q 



r 1- 



a 



(54) 



The result does not depend very much on the free parameter n. Roughly speaking, for the fiducial example discussed 
previously, the energy density of the particles becomes dominant 5 e-foldings before the exit of the horizon. Since 
the mode considered leaves the horizon 53 e-foldings before the end of inflation, we conclude that a model with ~ 60 
e-foldings does not suffer from the back-reaction problem. 

In the following, we turn to the calculation of the matter power spectrum taking into account the constraint derived 
above, namely a > 2. 
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FIG. 3: The function n(k) for fcb = 312, n = 1, a — 2 (continuous line), a = 5 (dotted line) and a = 10 (dashed line). 

2. Calculation of the power spectrum 

The first step is to calculate the two-point correlation function of the Bardeen potential. Most of the calculation has 
already been done in the previous subsection, see Eq. (36), but what matters now is the coefficient of proportionality 
which was not determined previously. One finds 



f 2 Q f°° r\k sin kr 

(* 2 (n,fc b )l*(»?.x)$(t ?) x + r)|* 2 (n,fc b )) = ^^ : / ™ S -^k 3 \f k \ 2 [l + 2nh(k)]. 



% 4tt Jo 



(55) 



The link between the power spectrum and the classical Fourier component of the Bardeen potential $(77, k), with 
$(77, x) = l/(2vr) 3 / 2 / dk$(?7,k)e 4kx , is obtained if one uses an ergodic hypothesis and identifies the "quantum" 
two-point correlation function with the spatial average ($(77, x)$(?7, x + r))v- In this case, one finds 



($(?7,x)$(77,x + r))v = 



dk sin kr , 



1 



; /r -fc3| $(??)k) |2^_ k 3 mvM 2 = ^_ Ask n s ~l [l + 2n - hm (56) 



Then, the matter power spectrum can be directly derived since the density contrast is linked to the Bardeen potential 
by the perturbed Einstein equations. As we mentioned above, we take n s = 1 for simplicity and get from the Poisson's 
equation 



l^7,k)| 2 = 



4 / ^phys 



|$fa,k)| 2 =H*(»7,k)| s 



7T £ 2 

— ^-A s k phys [l + 2nh(k)}, 



(57) 



where ^4 S is the constant fixed by the COBE normalization, see Eq. (39). The quantity (5(ry,k) in the last equations 
is dimensionless. The dimension-full (physical) Fourier component of the square of the density contrast is just Oq 
times the previous one, ao being the value of the scale factor today. Defining, as usual, the spectrum P(k) by 
P(k) = |<5 P hys(?7, k)| 2 /ciQ and choosing the normalization of the scale factor as ao = 2£u where £h is the Hubble radius 
today, one finds 



P(k) 



16tt 



'rms-PS 



1 



rp2 



W 9 2 



I 



— +2n^ 2) (fc b ) 

D7T 



1 + 2nh(k) 



phys 



(58) 
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This equation gives the initial matter power spectrum. In order to obtain the matter power spectrum today, we need 
to take into account the transfer function T(k) which describes the evolution of the Fourier modes inside the horizon. 
In that case, one has 



m=T\k)— A 

' ' n 



2 16-7T Q rms _ PS 1 



T 2 



1 + 2nh(k) 



^phys 



The sCDM transfer function is given approximatively by the following numerical fit [40] 



T(k) 



ln(l + 2.34g) 
2Mq 



1 + 3.89q + (16. lg) 2 + (5.46<?) 3 + (6.71g) 4 



-1/4 



, q = fc/p^Mpc- 1 ] 



where T is the so-called shape parameter, which can be written as [40] 



(59) 



(60) 



(61) 



where fio is the total energy density to critical energy density ratio and f2b represents the baryon contribution. More 
explicitly, the definitions used in this article are f2o = ^A + &m = + f^cdm + ^b- We have now normalized the 
matter power spectrum to COBE. It is important to realize that the above procedure only works for the sCDM model 
since we have used the Sachs- Wolfe equation (23). The sCDM matter power spectrum is depicted in Fig. 4. The 
measured power spectrum of the IRAS Point Source Catalogue Redshift Survey (PSCz) [41] has also been displayed 
for comparison. One notices that the effect of the step in h(k) is to reduce the power at small scales which precisely 



icr 




0.001 0.010 0.100 1.000 10.000 

k (h Mpc~ ! ) 



FIG. 4: Matter power spectrum normalized to COBE for different numbers n of quanta in the initial state. The cosmological 
parameters are those corresponding to the sCDM model, namely, h — 0.65, £7a = 0, fib = 0.05, f2 c dm = 0.95 and n s — 1. 
The parameters describing the non vacuum state are fc^ hys = 0.052/iMpc~ 1 and a = 2.5. The data points represent the power 
spectrum measured by the PSCz survey [42] . 

improves the agreement between the theoretical curves for n ^ and the data. Let us remind at this point that 
the shape of the function h(k) has not been designed for this purpose and comes from different (theoretical) reasons. 
Therefore, it is quite interesting to see that the power spectrum obtained from our ansatz fits reasonably well the 
data. This plot also confirms the result of Ref. [8], namely that the number of quanta must be such that (for sCDM) 
1 < n < 10, i.e., cannot be too large. 



14 



Another more accurate test to check the consistency of the model with the observations is to compute the rms mass 
fluctuation at a scale of r = 8/i _1 Mpc. Its definition for a top hat window function reads 



/ SM 



(r = S/i^Mpc) = 



1 



(4^/3)2 



with 



W(k) = Airrl 



sin kro cos krg 



(kr ) 3 (fcr ) 5 



(62) 



(63) 



For n = (i.e. standard sCDM model) with the following choice for the cosmological parameters, h = 0.65, O = 1, 
fib = 0.05, a numerical calculation of the previous integral gives as ~ 1.67, in agreement with previous estimates, see 
Fig. 15 of Ref. [43]. For h = 0.5, the result becomes <r% ~ 1.28 which is also consistent with Ref. [43]. Let us now 
calculate the rms mass fluctuation for n^O. For a = 2, fc£ hys = 0.052ft, Mpc -1 one finds, for our fiducial choice of 
the cosmological parameters, as — 0.99,0.79,0.68 for n = 1,2,3 respectively. In the same conditions, but for a = 5, 
one obtains as — 0.98, 0.77, 0.66 which illustrates explicitly the fact that er 8 is not sensitive to the parameter a. These 
numbers should now be compared with those inferred from cluster abundance constraints [44] which give 



erf 1 = (0.5 ±0.1)O" 



(64) 



where 7 ~ 0.5. We recover the well-known conclusion that the n = sCDM model is ruled out because the CMB and 
clusters normalization are not consistent with each other (i.e. the difference is greater than 5a). On the other hand, 
we see that putting a few quanta improves the situation and that the n = 2 model becomes compatible with the data 
at the 3a level whereas the model n = 3 gives the correct value at less than 2a. Fine-tuning the other cosmological 
parameters would allows us to achieve an even better agreement. Therefore, as announced, there exists a region in 
the space of parameters where the model is in agreement with the presently available data. 

Now, we would like also to test the predictions of the model in the case where a cosmological constant is present. A 
first problem is that the value of the coefficient ^4 S is no longer the same. The reason is that the integrated Sachs- Wolfe 
effect now plays a role at large scales and modifies the relation (24). This changes the COBE normalization and we 
have C^^Oa 7^ 0) 7^ C^(fl\ = 0). As a consequence, the constant ^4 S in Eq. (26) is no longer given by Eq. (39) and 
has to be evaluated numerically for each value of the cosmological constant. We can parameterize this dependence by 
introducing a coefficient B s (fl\) such that ^4 s (Oa 7^ 0) = B s (Oa)A s (Oa = 0). Obviously, one has -B s (Oa = 0) = 1. 
A second problem is that the transfer function is also modified by the presence of a cosmological constant. However, 
the change can be easily parameterized by using a new numerical fit to the transfer function. Finally, the power 
spectrum can be written as 



P(k) = T\k) 



2 /,^3 2 (^o) D ^ ^ 16tt <2? ms _p S 1 



5 2 (a, 



-B S (Q A ) 



5#o 4 



T 2 



1 

6tt 



+ 2nD^\k h ) 



1 + 2nh(k) 



phys 



(65) 



where the function g(fl) takes into account the modification induced in the transfer function by the presence of a 
cosmological constant. Its expression can be written as [45]: 



50 
~2~ 



4/7 - fl A + [ 1 



1 + M 

70 J 



Using the previous expression, we deduce that 



er|(O , O c dm, Ob, Oa, n) 



.9 2 (Oo) 
.9 2 (O m ) 



£? s (Oa)<7 2 (O , O cdm , O b , Oa = 0, n), 



(66) 



(67) 



where the last term, cr|(Oo, O c dnn Ob, Oa = 0,n) does not depend explicitly on O ct j m since this one does not appear 
in the shape parameter and therefore neither in the transfer function. Therefore, this term is equal to the value 
computed previously for the sCDM model. From this last equation one can easily deduce the constant £? s (Oa), and 
we finally obtain 



cr|(O , O c dm, Ob, Oa, n) 



g|(Oo, Ocdm, Ob, Oa, n = 0) 
erf (O , O c dm, Ob, Oa = 0, n = 0) 



oJ(O , O cdm , O b , Oa = 0, n) 



(68) 



The quantity cr|(O , O c d m , Ob, Oa, n — 0) is known in the literature since it is nothing but the value for the ACDM 
model. For our choice of the cosmological parameters, one has, for the ACDM model with Oa = 0.7, as — 1 in 
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agreement with Fig. 16 of Ref. [43]. The quantity c7g(f2o> ^cdmi ^t>; ^A = 0,n = 0) is the value corresponding to 
sCDM. Finally, the quantity <7g(Qo> ^mi ^t>) = 0, n) is known from the previous calculations. Hence, for the values 
of the parameters considered above, one finds erg ~ 0.59, 0.47, 0.40 for n = 1, 2, 3 respectively. This has to be compared 
with tTg lu which is now equal to CTg lu ~ 0.91 [46, 47, 48]. In this case, we see that already the first value n = 1 gives a 
too small contribution. This had already been noticed in Ref. [8] for the CMB anisotropy. In this case the presence 
of the cosmological constant increases the height of the first acoustic peak which is also the effect of adding more 
and more quanta. This can result in an acoustic peak which is too high. The cure is obvious: one has to decrease 
the value of the cosmological constant as also noticed in Ref. [7]. For instance, for f^A = 0.5, the value of as for the 
ACDM model becomes a$ ~ 1.25, see Fig. 16 of Ref. [43]. This gives for our model cr 8 ~ 0.74, 0.59, 0.51 for n = 1, 2, 3 
which goes in the right direction. The model with n = 1 is now compatible at the 1.5cr level. One can decrease more 
the value of Oa in order to obtain a better agreement. These numbers are in full agreement with those proposed in 
Ref. [7] where the range for the cosmological constant in a BSI model with p > 1 was found to be 0.2 < Cl\ < 0.5. 
Again, this is not surprising since we essentially deal with (almost) the same two-point correlation function. 



IV. FOUR-POINT CORRELATION FUNCTION FOR NON- VACUUM INITIAL STATES 

A. General expressions 



In this section we proceed with the calculation of the four-point correlation function. We first perform the calculation 
for the state which we then generalize for the state l^)- The first step is to establish the expression of all the 

combinations of four creation and/or annihilation operators taken in the state \^\(a,n)). One finds 



(*i(CT,n)|ckCpctcl|*i(CT,n)) = 



(*i(tr, n)|4ct,c q c s |*i(CT, n)) 



(*i(cr, n)|cj c c P 4c s |1'i( ( 7, n)) 



n)|c k c p c q 4|*i(cr, n)) 



6(p - s)8(k - q) + 6(p - q)6(k - s) 

1 + nd(s e V) + n8(q G V) + n 2 S(s G V)5(q G V) 
n(n + l)S(s G V)5(q - s)S(p - s)<5(k - q) , 



= n 2 5(s G 2?)<5(q G T>) 



S(p - s)<5(k - q) + 5{p - q)5(k - s) 



-n(n + l)6(s G V)5(q - s)6(p - s)<5(k - q) , 

n6(s G V)6(p - q)6(k - s) + n 2 6{s G V)6(q G V)6(p - q)5(k - s) 

+n 2 S{s G V)6(p G V)6(q - s)5(k - p) 
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+nS(q G V)8(p - q)6(k - s) + n6(p G P)<5(q - s)5(k - p) 



+n 2 S(s G 2?)(5(p G D) 



J(q - s)<5(k - p) + 6(p - q)8(k - s) 



-n(n + l)6(s G V)8(q - s)6(p - q)8(k - s) . 



(69) 
(70) 

(71) 



(72) 



(73) 



(74) 



16 



Then we can use the previous equations to calculate the expectation value of four coefficients ag m in the state 
\^i(a,n)}. Using Eq. (28) which links the operators ag m to the operators c k , one obtains 

£ 4 

{^i(a,n)\ae imi ae 2m2 ae 3m3 ae 4m4 \^i((7,n)) = 
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Let us notice the following technical trick. Originally, in front of the first squared bracket in Eq. (75) appears a term 
of the form (— l^i+^gW^i+^+^+M/^ Using the fact that the presence of the Kronecker symbols implies that the 
corresponding expression inside the squared bracket is non vanishing only if £\ = £3 and £ 2 = £4, the previous term 
can be rewritten as e 2l1T ( f ^+^) = 1. This explains why it does not appear explicitly in Eq. (75). Similar manipulations 
can be performed for the terms in front of the following two squared brackets. Let us also notice that Eq. (75) includes 
a complex exponential factor, namely exp[i7r(^i + £ 2 + £3 + &i)/2]. However, by inspection of the properties of the 
quantity W™ 1 ^ 2 ^ 3 ™ 4 , see its definition in terms of Clebsh-Gordan coefficients in the Appendix, it is possible to show 
that this term is in fact real. Indeed since the Clebsh-Gordan coefficients are non-vanishing only if £\ + £2 + L = 2p 
and £3 + £4 + L = 2q where p, q are integers, we have £\ + £2 + £3 + £4 = 2(p + q — L), i.e., an even number. Therefore, 
the complex exponential factor in Eq. (75) is real and equal to either 1 or —1. These technical considerations will be 
employed in the formulas below. We also see that we now have (—1)^1+^2+^3+^4 _ \ Finally, Eq. (75) can be cast 
into a more compact form 

(^i{(7,n)\a eimi a e2m2 a e3m3 a e4m4 \^ 1 (a,n)} = -j^-l {-l) mi+m2 C^c£ > 6e 1 t 3 St 2 t i Sm 1 -m 3 tim 2 -m 4 
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The first part of this equation has the same structure as the corresponding well-known equation for the vacuum state. 
It is sufficient to replace Ct with in the latter to obtain the first part of Eq. (78). However, in addition, there 

is a non trivial term proportional to the coefficient E^^^ which cannot be guessed a priori. Obviously, for n = 0, 
one recovers the standard result. 

The calculation of the four-point correlation function in the state ^2) is a bit more involved. Using the fact that 
the operators ag m do not act on ct, one finds that the general expression is given by 

1*00 

|* 2 (n,fc b ))= / dCT 3 2 (CT;fc b )(*i(CT,n)|a £l |*i(a,n)>. (79) 
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The integration of terms of the type Ce 1 Ci 2 , C^D^J and E^\ 2t ei is easy and proceeds as before. The most difficult 

part is the integration of terms of the type We find that, in the state ^2}, the four-point correlation 

function is given by 



(^2(n,k h )\ai imi ai 2m2 ae 3m3 ae 4m4 \^2(n,k h )) = 
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with 
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We now see clearly the complication brought into the problem by the term F^ ' t . This term prevents us to reduce 
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the terms within the squared brackets to the natural form cf^C^ because F^\ ^ D 
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Calculation of the excess kurtosis 



We are now in a position to calculate the excess kurtosis. In the previous section, we have established the expression 
of the four-point correlation functions for the operator ag m . In order to establish an analytical formula for the CMB 
excess kurtosis, one just needs to use the equation linking ai m and 5T/T and to play with the properties of the 
spherical harmonics. Explicitly, the excess kurtosis is defined as 

JC = fi4 - 3/4 > (83) 

where the second moment has already been introduced and where the fourth moment, /U4, of the distribution is defined 
as 
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fXi = (K) with K 



ST 
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(84) 



An important shortcoming of the previous definition is that the value of JC depends on the normalization. It is much 
more convenient to work with a normalized (dimcnsionless) quantity. Therefore, we also define the normalized excess 
kurtosis as 



M2 



(85) 



which is the one more commonly used in the literature. In what follows we work with either JC or Q parameters. 
Thus, Eqs. (33), (80) and (83) imply 



JC = 



'pi 



1% \ (4tt) : 
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(86) 



Let us first concentrate on the first term in the above equation. The terms (2£i + 1)(2* 2 + 1) and l/(47r) 2 originate 
from the addition theorem of spherical harmonics 



47T 



(87) 
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(2) (2) (2) 

The factor 3 comes from the definition of /C, see Eq. (83). The fact that F£ £ ^ D\ D\' prevents this first term from 
vanishing. This is consistent with the previous considerations, as we have seen that in the absence of this condition 
the structure of the four-point correlation function would be similar to the one in the vacuum state, up to the term 

( 2) 

proportional to ' ( ( ti of course. Let us now treat in more detail the second term in Eq. (86). Using again the 
addition theorem of spherical harmonics and the expression of a Lcgendre polynomial in terms of a spherical harmonic, 
Y m = v /(2£+l)/(4^)P £ (cos0), we can perform the sum over the indices m,'s and express the corresponding factor 
in terms of the coefficient Hg^i t , and therefore in terms of Clebsh-Gordan coefficients, see the Appendix. 

After some lengthy but straightforward algebra, one finds that the excess kurtosis in our class of models is finally 
given by 
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(88) 



Let us emphasize that Eq. (88) is the general expression for the excess kurtosis for any non-vacuum state, since the 
only information we have used about the function h(k) is that it is always positive, it vanishes at infinity, and it is 
a monotonically decreasing function of k. This expression is just a pure number, and it is our main result. In the 
following, as we did in previous sections, we will choose an adequate ansatz for h(k), namely that one from Eq. (43), 
and compute the excess kurtosis JC, as well as the normalized excess kurtosis Q defined by Eq. (85). We will then 
compare the calculated value for Q to the one quantified by the cosmic variance. 

In an analogous way as for the definition of the second moment \i 2 given in Eq. (42), and for future convenience, 
we can express the excess kurtosis in terms of its "multipole moments" K-i l i 2 g 3 i i , as 



I1I2I3U 

Then, from Eqs. (34), (81) and (82), it is easy to establish that the moments K^t 
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can be put under the form 
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(90) 



The last step consists in normalizing the spectrum. For that, we use the value of A s determined previously (in the 
sCDM case with non- vanishing quanta n in the vacuum state). We obtain 
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In particular, we have the following expression for /Q 
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This expression for the multipole moments will be employed in the next section to estimate the overall amplitude of 
the non-Gaussian signal from non- vacuum states in a semi-analytical manner. 

Let us end this section by signaling that the explicit expression for the normalized excess kurtosis parameter can be 
easily derived from the above formulas. Then, since this derivation is not especially illuminating, we prefer to jump 
directly to the numerical evaluation. 



V. RESULTS 



Having established the formal expression of the excess kurtosis, we now turn to the question of its numerical 
evaluation. As we are going to see, it turns out that it is not possible to calculate everything analytically for our 
specific ansatz. Therefore, we start with giving a detailed order-of-magnitude estimate of the excess kurtosis employing 
the model described above. After this, in the second subsection, we present a full numerical evaluation of Q. The 
analytical estimate will help us in roughly understanding the full numerical results. Finally, a comparison with the 
cosmic variance of the excess kurtosis will tell us about the feasibility of detecting this non-Gaussian signal. Readers 
interested in just the final result can skip this first subsection and continue reading in section V.B. 



A. Approximate analysis 



In what follows, we estimate the order of magnitude of the expected excess kurtosis. We first study the function 
-D^(cr), defined in Section III as 



Df \a)^ ^ J, 2 +1/2 (fc)fc" s - 3 dfc . 
Jo 



(93) 



The plot of this function is represented in Fig. 5. We can easily understand the qualitative behavior of this function. 
For small values of the argument, we take the first term of the Taylor expansion of the Bessel function and perform 
the integration exactly. The result reads 
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7 M +»S-1 (T<< 1 



2 2 ^+ 1 (2^ + n s - l)r 2 (^ + 3/2) 
On the other hand, for large values of the argument, using Eq. (6.574.2) of Ref. [39], we obtain 

r[3-ns]r[*+(ng-l)/2] 



23-n sr 2[2 - n s /2]T[£ + (5 - n s )/2] : 



a > 1 . 



(94) 
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For n s = 1, the above amounts to Df\a) ~ l/[n£(£ + 1)] which, for £ = 2, gives D^'(a) ~ 1/(6tt) ~ 0.053 in 

- (2) 

agreement with Fig. 5. As a next step we want to understand the qualitative behavior of D e as a function of fcb and 
£. Its definition, given in Section III, reads 
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where we have used the fact that the behavior of h(k), especially for large value of the parameter a, is very similar 
to a Heavisidc function. In practice, for large fcb and for the range of angular frequencies we are interested in, we will 

— (2) — (2) 

have fcb ^> £■ In this case, D\ ~ 1/[tt£(£ + 1)]. Let us now try to evaluate F£ £ defined by 
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Using again the fact that h(k) behaves like a step function, we find 
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FIG. 5: The function Df\<j) versus a for I = 2. 



The previous equations allow us to estimate the first term in Eq. (88). We find 

F?X-T)fDf}^. (99) 

This conclusion rests on the approximation that the function h(k) behaves like a step function. In reality, this is of 
course not exactly the case. It would have been nice to find the order of magnitude of the correction in terms of 
the parameter a controlling the sharpness of the function h(k). Unfortunately, the complexity of the equations has 
prevented us from deriving such an estimate. This notwithstanding, we can be sure that the contribution of this first 
term to Eq. (88) will not be higher than that of the second term, see below. 

The second term which participates to the excess kurtosis is given by the term proportional to the coefficient 
Ejfy 2i i4 - Let us now study this coefficient in more detail. Its expression, using the approximate behavior of h(k), 
can be written as 

/■+oo 

E?Xe 3 u = / 4 + i/2(^)4ti/ 2 W4+i/2(^)4 + i/2(*(^ 2 "" 8 ^ (ioo) 
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* l " ^ 1 + l/2(fc)^ 2 + l/2(fc)^ + l/2(fc)^4 + l/2(fc)fc 2nS_8 dA : . (101) 
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We have not been able to find a compact expression for this coefficient. However, it is easy to follow qualitatively 

its behavior. Its numerical value is maximum when the indices £j's are all equal. This is because this configuration 

maximizes the overlap of the first peak of the Bessel functions. Also when t\ = £2 = £3 = f-4 increases, the numerical 
— (2) 

value of Ef' iei decreases, because the coincidence between the four (first) peaks of the Bessel functions occurs 
at higher k, resulting in a factor k 2ns ~ 8 much smaller. From these considerations, one can infer that the largest 
coefficient is £2222 • Practically, this quantity does not depend on fcb because we always choose kb 3> 5/2 so that the 
upper bound of the integral can be considered to be the infinity. Numerically, one obtains 

E { 22 \ 2 = 0.000094 ~ 10~ 4 . (102) 

This means that for all £j's we have £^] 2 ^ 4 < 10~ 4 . 
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These semi-analytical considerations allow us to derive a rough estimate of the excess kurtosis. In the context of 
our approximation, this one does not depend on fc b . Let us first write Eq. (88) as follows: 



/C = ^(/C 1+ /C 2 ) 
and consider the first term. One has, as already mentioned previously, 
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Let us now study the second term. Since we have shown that the contribution of -E 2 222 dominates the sum, we can 
very roughly estimate the excess kurtosis by retaining only this term. We have 
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where we have used the fact that W2 — 1 for the COBE window function, see below. The sum of the Clebsh-Gordan 
coefficients can easily be computed by noticing that 
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the other coefficients being zero. The numerical value of the sum is 0.0857 and we therefore reach the following result 
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Now, to have an order of magnitude of the contribution to the kurtosis from the second term, JC2, what remains to 
be done is to take into account the normalization. Using the expression of A s derived previously into Eq. (103) we 
finally find 



JC ~ -10~ 2 7T- n(n + ^ " ^rms-PS 



(2n+l) 2 
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(109) 



Let us again stress that there is a non-trivial guess in this calculation that can only be justified by the full numerical 
calculation. The contribution of JC2 is small, JC2 <C 1. Since we have taken K.\ ~ 0, this means that, in fact, we 
have assumed /C2 ^> K,\. This has clearly be done without a proof because we were not able to derive an order of 
magnitude of K\ as a function of the sharpness parameter a. Nevertheless, we will see that the previous estimate is 
quite good. 

It is also very convenient to work with the normalization-independent quantity Q. Recalling that with our approx- 
imations we have the second moment 
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it is easy to show that Q (= Qi + Q 2 — Q2) can be written as 
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In the above equation we have used the COBE-DMR window function Wt ~ exp \—^£(£ + 1)(3.2°) 2 ] , where 3.2° is the 
dispersion of the antenna-beam profile measuring the angular response of the detector and we have chosen £ max = 20; 
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we find that the sum in the expression for /x| is equal to ~ 3.73. As an example let us take n = 2. Then we find 
Q ~ —4.09 x 10~ 4 . Of course this number should be considered only as an order-of-magnitude estimate. It can be 
easily improved if we add more terms in the calculation of the sum in the term K.2- At this point, one should make 
clear that taking into account more terms does not mean that we consider the next-to-leading order of a consistent 
expansion. In the present context there is no small parameter to expand in. Therefore, the choice of the extra terms 
that we include in the sum is a bit arbitrary. However, this is not a serious problem since we know that the term 

(2) (2) 

Ef (( ( gives a contribution greater than the contribution coming from E e , e , e , £ , provided that £j -C l\, i = 1 • • • 4. 
In other words there is no mean to give a precise rule with regards to the terms that should be kept or not kept but 
there is clearly a general tendency which renders the improvement of the approximation possible (otherwise the sum 
would not be convergent). A good strategy to calculate the successive corrections is to proceed as follows. Suppose 
that we would like to calculate the sum K-2 up to £ max . A consistent requirement is that one should take into account 

(2) (2) 

all the terms E\' ttt such that £i < £ max , i = 1, • • • 4. For example if £ max — 2, then one should only include ££222 as 
we did in our "leading order" calculation. The next step is to consider the case £ max — 3, i.e., we choose to take into 

— ( 2) — ( 2) — ( 2) ~ (2, ) 

account the terms £2233: -^33331 -^2333 an d £2223 (°f course we should also include the terms obtained by permutations 
of the indices £,). In fact, it is easy to see that the last two terms do not give an extra contribution because the sum 

^ —(2) — C2) 

Y^i^-i is n °t even. For the two other terms, we find £2233 — 0.000041 and -E3333 — 0.000021. After a lengthy but 
straightforward calculation, one can show that this gives an additional contribution equal to —8.51 x 10~ 6 . Therefore, 
a more accurate estimate of the normalized excess kurtosis is 

Q ~ -4.17 x 10~ 4 . (112) 

We now turn to a full numerical computation of this quantity. 

B. Full numerical results 

In the previous section we have developed a rough idea of the amplitude of Q. This gives us an estimate of the 
excess kurtosis parameter. However, due to the various approximations, a full numerical resolution of Eq. (85) proves 
essential. 

By resorting to Eq. (89) and using again the COBE-DMR window function, we can compute the value of Q, valid 
on large angular scales. We plot the results in Fig. 6, where we show the normalized excess kurtosis Q for some 
particular values of the free parameters. As we see from it, Q ~ —4.24 x 10~ 4 is an asymptotic value, provided 
we concentrate on the middle and big values of the built-in scale kb- This value almost exactly corresponds to the 
numerical estimate derived previously. In other words, the fact that the numerical estimate does not depend on fcb is 
confirmed by the plot, except for small values of the wave-numbers. In fact, this shows that the quantity Q does not 
depend very much on the free parameters. We have already established this property for a but this is also true for n 
since, using the analytical estimate, we find, for n — 1, Q ~ —3.77 x 10~ 4 and for n — > 00, Q ~ —4.25 x 10~ 4 . Since 
we know that this result does not depend on the details of the weight function h(k), we conclude that the asymptotic 
value obtained above is a generic value, at least for large values of kb- In particular, this is true for kb ~ 300 which 
corresponds to the built-in scale located roughly at the privileged scale in the matter power spectrum selected by the 
redshift surveys of Ref. [27]. Another important remark is that the excess kurtosis is found to be negative. 

Let us now try to understand qualitatively the shape of the plot Q vs. kb- The state that we consider, 1*2), is a 
quantum superposition of states l^i), each one of these containing rife quanta for all the scales k up to a given, fixed 
scale o [cf. Eq. (17)]. The "weight" given to each state |\f f i) is described by the function g(k; kb) that depends on the 
privileged scale kb- As already mentioned, for our ansatz of Eq. (43) we have g(k; kb) — (a/2fc) 1 / 2 cosh -1 [a ln(fc/feb)], 
which, as a function of k, is roughly "peaked" at kb- Then, in effect, we may approximately write 

|* 2 (n,fcb)>^|*i(n,fcb)>= (£) l nfe > I°p>- ( 113 ) 

kev(k h ) P ^D(k b ) 

Thus, we see that a small fcb will reduce the range of scales k included in the domain T>{kb), and therefore also reduce 
the effective available number of quanta (of energy k). Given that we employed the Sachs- Wolfe formula, the excess 
kurtosis that we computed is only the "trace" of the non-Gaussian signal characterized by kb left at large scales [this 
is also why there is no contradiction in using Eq. (24) while kb can be large]. Let us consider two scales, say, kbi and 
kb2 such that kb2 > febi ^> 1 (recall that we are taking 770 — i]i ss = 1)- It is clear that passing from kbi to fcb2 will not 
change the structure of the state ^2) at small scales. It will just enlarge the domain where there are quanta, leaving 
unmodified the large scale part. Therefore, as the excess kurtosis is essentially given by the large scale part of [^2), it 
must be independent of fcb provided that kb is large, exactly as we find. On the other hand, if kb is small, a change in 
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FIG. 6: The normalized excess kurtosis parameter Q in terms of the privileged (comoving) wave- number fcb for a particular 
representative set of parameters: a = 5 for the sharpness of the weight function signaling the privileged scale, and n = 2 quanta 
in the non-vacuum initial state for the cosmological perturbations. 



this scale affects the structure of the state at scales that are relevant for the Sachs- Wolfe effect. As we saw previously, 
in this regime, the number of quanta corresponding to large scales decreases as goes to zero. The result is that 
the excess kurtosis, whose value is directly dependent on the number of quanta, should diminish proportionally as 
goes to zero, and this is in fact what we see in Fig. 6. 

These considerations give an intuitive understanding of the main result of this article. We will see in the next section 
that the theoretical uncertainties on Q (as given by the cosmic variance) are roughly equal to Q cv ~ 1 implying that 
the non-Gaussian feature studied previously would be undetectable. 



C. Comparison with the cosmic variance 

The cosmic variance quantifies the theoretical error coming from the fact that, in cosmology, observers have only 
access to one realization of the ST/T stochastic process whereas theoretical predictions are expressed through ensemble 
averages. To compute the cosmic variance of a given quantity, say e, one should proceed as follows [6, 52]. Firstly, 
one has to introduce a class of unbiased estimators E of the quantity e, i.e., (E) = e, where in the present context 
the average symbol means a quantum average in the considered state. Secondly, one should compute the variance of 
these estimators and find the smallest one under the constraint that the estimators are unbiased. The estimator which 
possesses the minimal variance is the best unbiased estimator of the quantity e. We denote it as E^ est . Finally, one 
should compute the variance of the best estimator up which is, by definition, the cosmic variance. If up = 0, 
then each realization gives e and from one realization we can measure the quantity we are interested in. If <7g ^ 0, 
which is obviously the usual case, one can attach a theoretical error to the quantity we seek. Let e rca i be (the numerical 
value of) one realization of the corresponding stochastic process, then we can say that the quantity e is found to be 
e rca i ± erg So far, this strategy has been successfully applied to quantities related to the two-point correlation 
function [6] and to the three-point correlation function [52], in the vacuum state. 

The normalized excess kurtosis is defined in Eq. (85). We see that, in the present context, we face two important 
complications. The first one is that one has to perform the minimization in a non-vacuum state. Let us notice 
in passing that this implies that the cosmic variance of the multipole moments is probably not given by the usual 
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expression if the quantum state is no longer the vacuum. The second complication is that, in order to determine the 
best estimator of the normalized excess kurtosis, one has to deal with the ratio of two stochastic processes. Suppose 
that we want to find the best estimator of the quantity e = ei/e 2 knowing the best estimators of the quantities ei 
and e 2 , i?best(ei) and E hcst (e 2 ), respectively. The problem is that E hest (e-L / e 2 ) ^ E^^e-i) / E hest (e 2 ) . In this case, 
the calculation of the best estimator becomes much more complicated. 

Therefore, the full calculation of the best estimator of the normalized excess kurtosis is a project beyond the scope of 
the present work and we will instead limit ourselves to an order-of-magnitude estimate of the theoretical uncertainties. 
This will yield a roughly correct estimate of the uncertainties without the complications of much more cumbersome 
analysis. 

To avoid the first complication, we estimate the variance of the excess kurtosis as if it were issued from a Gaussian 
process, i.e., as if the quantum state were the vacuum state. To deal with the second complication, we use the 
following procedure. The excess kurtosis JC (or its normalized version Q) is just related to the fourth moment /X4 
of the distribution from which we substract the Gaussian part ^ auss ^ = 3^i 2 . We will only take into account the 
contribution to the cosmic variance related to the fourth moment /14. This last one is determined in the standard 
manner: the quantity K, defined in Eq. (84), is an unbiased estimator of /i 4 ; let <t cv be its variance (in order 
to compute this variance one needs to calculate the eight-point correlation functions of the relevant creation and 
annihilation operators in the vacuum state-according to our trick to avoid the first complication, see above-) and let 
^roai _ ^^reai one reauza tion of the stochastic process K; as it is c cv what attaches theoretical error bars to the 
actual value for the mean kurtosis, we can heuristically express the effect of <r cv on ^4 as follows: /14 ~ /i™ 3,1 ± <r cv , at 
one sigma level. Having specified the cosmic variance of the kurtosis, we now need to relate it to the cosmic variance 
of the excess kurtosis. The effect of the cosmic variance is that, instead of finding the value Q = for a Gaussian 
process in the vacuum state, we typically obtain a value shifted by ± Q cv which can be estimated as 

^ ^4° al ~ 3^2 ^ M4 ~ 1 a cv : a cv (-\-\a\ 

Jcv ~ A ~ Mi Mi ~ Ml' { ' 

where in the last equality, and as we mentioned above, we used the fact that the process is Gaussian. In other words, 
Eq. (114) shows that Q cv is the normalized excess kurtosis parameter (assuming Gaussian statistics) purely due to 
the cosmic variance. Q cv is in general non-zero [49] and its magnitude increases with the theoretical uncertainty 
ct cv . This gives a fundamental threshold that must be overcome by any measurable excess kurtosis parameter. The 
expression for Q cv is straightforward to obtain, although after a somewhat long algebra. It was computed in Ref. [34] 
and gives 



E^ + i)c £ wf] 2 

E*(2* + i)cW J ■ ( j 

As is the case for Q, one of the advantages of the previous expression for Q cv is that it is transparent to any particular 
normalization of the spectrum. We have computed the value for Q cv for the COBE-DMR window function and found 
a value of order one for values of the scalar spectral index close to one. For illustrative purposes, we show in Fig. 7 
the variation with spectral index of Q cv when the term with pre-factor 24 is absent [just a 5% off of the full result]. 
Note that in an analogous calculation in Ref. [34] the quadrupole was subtracted from the sum. The fact of including 
it now only increases slightly the final value for Q cv , showing the big contribution of the low order multipoles to the 
cosmic variance, as already noted in that paper. 



VI. CONCLUSIONS 



In this article, we have presented evidence that Gaussianity is a robust property of inflation. We have seen that a 
departure from the standard vacuum initial conditions for the cosmological perturbations leads generically to a clear 
non-Gaussian signature, viz. the excess kurtosis of the CMB temperature anisotropies. The signal-to-noise ratio for 
the dimensionless excess kurtosis parameter is found to be 

~4xl(T 4 , (116) 

and so the signal lies well below the cosmic variance and far away from experimental detection. We have found that 
this value is quite independent of the free parameters of the model. We have also shown that the excess kurtosis 
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FIG. 7: The normalized excess kurtosis parameter Q cv due to the cosmic variance, in the mildly non-Gaussian approximation 
and as function of the scalar spectral index n s . We have checked that, of the two terms in Eq. (115), it is in fact the first 
one (with pre-factor 72) the one that dominates the full result for Q cv . It is this first term what we plot in the figure for the 
COBE-DMR window function, including the quadrupole and for i? m ax = 20. The addition of the term with pre-factor 24 raises 
the points no more than a 5% for all values of ns . 



is generically negative. The only possible loophole in the argument presented above, and that we have discussed at 
length, is the uncertainty related to back-reaction. This issue is generically important for the inflationary phase, and 
one could well conceive that it could modify the evolution of the background in such a way as to increase the ratio 
S/N. However, we do not think at present one such positive conspiracy would take place, and therefore primordial 
Gaussianity keeps on being a generic property of single field inflationary models. 

Finally, a comment is in order on the trans-Planckian problem of inflation. As already mentioned, it has been 
recently suggested in Ref. [17, 18] that trans-Planckian physics could be mimicked by a non-vacuum state (rather 
than by a change in the dispersion relation) implying that non-Gaussianity would be a possible observable signature 
of the physics on lengths much smaller than the Planck length. Unfortunately, what has been shown in the present 
paper for the CMB excess kurtosis leads us to conclude with a no-go conjecture: these aforementioned signatures of 
trans-Planckian physics will most probably be astrophysically unobservable. 
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Appendix 



We review here some definitions and properties of quantities involving Wigner 3-j symbols which are useful for the 
main text. We define the integral of four spherical harmonics as 
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where, following the notation of Ref. [50, 51, 52], we wrote 



nZWTir = n Z\7r = J W W s ( e ) • ( 119 ) 

^t&tl™' 3 nas ^ ne following simple expression in terms of Wigner 3-j symbols [53]: 



^ mi m 2 m 3 _ Wl + 1)(2^2 + 1)(24 + 1) ( h *2 h\( h ?2 h , 

rt ^«3 - V 4^ ^ )\ mi m 2 m 3 J ' U j 

In the main text we employed the quantity H™^ 2 ^ 3 ™ 4 which is simply related to U^^lT*™* defincd above b y 

q,mim 2 m 3 m4 _ ( \£\+£ 2 \nDm 1 m 2 m 3 m l -i* _ ( \ti+l 2 'Qm 1 m 2 m 3 m i ;ini\ 
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The last equality holds since the 3-j symbols, as well as the Clebsh-Gordan coefficients, are all real. 



